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ABSTRACT 



The measurements described represent the first step in an experi- 
mental test of a classical impulse approximation for the capture process 



H + CH, 



H + CH„ . Since the predicted cross-sections in the 
2 3 



”18 2 

applicable energy range are on the order of 10 cm or less, the 
focussing properties of a thin lens beta ray spectrometer will be 
utilized to obtain maximum collection efficiency for the ions. 

The result of measuring the axial component of the magnetic field 
strength is given. A current-carrying wire was used to simulate the 
trajectories of the ions in the energy region 50-900 ev originating 

near the center of the magnet with a departure angle between 35 and 55° 
with respect to the magnet axis. Measurements of the trajectories as 



well as the variation of departure angle with magnetic field strength 



at a fixed energy, 



36 . a . 

T'tr , were carried out. 
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Scattering angle in laboratory frame, azimuthal angle used in 
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Introduction 



The data and results presented in this paper are part of an experi- 
mental project presently being conducted at the U, S. Naval Postgraduate 
School to verify a classical theory developed by Bates et, al. (1964a) 
for the capture of a light ion or atom from a target system by a fast 
projectile , 

The particular reaction selected to be experimentally investigated 
is the capture of a hydrogen atom (H) by a fast proton (h"*~) in its 
passage past a molecule of methane (CH^) „ This can be conveniently 
written as 

H + + HCH 3 - H + H + CH 3 . (1) 

The energy range over which a classical description of this reaction 
is valid is examined . The de Broglie wavelength and Heisenberg uncertainty 
criteria, plus the additional condition that the energy of motion of H + 
relative to H must be less than the binding energy of the combined 
product , limit the energy range to about 25-800 ev, and the scattering 

angle to the vicinity of 45 degrees. 

The standard procedure of measuring the number of particles entering 

a detector placed at a definite scattering angle will not be feasible 

- 18 2 

because of the expected low cross section (10 cm or less). In order 
to obtain maximum counting efficiency, a thin lens beta ray spectrometer 
was selected which permits all particles coming out in a cone at a given 
scattering angle to enter the detector. 

The axial component of the magnetic field was measured both along 
and off the axis with the results tabulated in Appendix 1. As was 



. 1 . 



expected, the field strength is linear with magnet current and is 
largest at the center,, In order to provide insight into the positioning 
of the interaction chamber, the detector, the various baffles, it was 
necessary to have &ome information concerning the trajectory of the 

particle,. This was done by simulating the particle path by a flexible 
current-carrying wire in the magnetic field , It is shown that the wire 
takes up a position coincident with the orbit of the charged particle 
provided that the current and tension of the wire are adjusted to 
correspond to the energy of the particle in question through the 
expression 

E^ 2 = 48 r (2) 

1 



where E is the energy of the particle (also that of the incoming 

“f- 

H ), M is the tension of the wire in mass units, and i is the wire 

current. Wire trajectories corresponding to particles having energies 

from 50~ 900 ev and leaving the magnet axis with angles from 35-55 degrees 

are examined and tabulated in Appendix 4. The variation of the particle 

departure angle with respect to the magnetic field strength at constant 
d9 



energy 



dE 



is also calculated from the experimental data. All of 



the trajectory work was done with the wire leaving the axis near the 
center of the magnet and focussing at a fixed distance approximately 
40 cm down axis. It was experimentally determined that optimum trajec- 
tories occurred between these two points. 

At the present time a theoretical analysis of the particle trajec- 
tories in the same energy and angular ranges is in progress. It is 
hoped that these numerical computations will substantiate the experimental 
data given in this paper. 
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2. Theory 



(a) Impulse Approximation 

Since a quantum mechanical treatment is not possible for any 
but the very simplest atomic interactions, the use of approximate solutions, 
particularly semic lass ica 1 , become important. 

The impulse approximation states that for high energies where polar- 
ization forces and chemical binding of the colliding particles may be 
ignored, the components of a complex system interact independently of 
each other. Furthermore, under certain conditions, classical theory 
may be used to examine the interaction and predict cross sections. 

We shall be concerned with the following ion-molecule rearrangement 
collision which can be handled experimentally 

H + + CH 4 - H 2 + + CH 3 (3) 

and more conveniently written as 

H + + HCH 3 h\ + CH 3 

where the underlined terms represent particles having considerable 
kinetic energy in the laboratory frame. 

Using Fig. 1 as portraying the capture mechanism, the reaction is 
briefly described in the following manner. H (a hydrogen atom) initially 
bound to CH^ (a methane molecule) suffers a binary collision with H + 

(the proton) and then suffers another binary collision with C (the carbon 
atom). Classical kinematics predicts that H + and the doubly scattered 
H will have a small relative velocity such that there is a large 
probability that the two particles will be bound and emerge as an 



particle . 
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Fig. 1 

Hydrogen Atom Capture Mechanism in Laboratory Frame 
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If the assumption is made that the nucleus of C is infinitely 



heavy, then in the laboratory frame 



ej 45° , 9^ ~ 45° , 9£ 90° 



and 



v ' =* v 1 v " 

1 2 2 




(4) 



It is important for both theoretical and experimental reasons to 
determine over what energy range this classical concept is valid. The 
requirements that must be satisfied are the usual de Broglie wavelength 
limitations and Heisenberg uncertainty relations for both binary collisions 
plus the additional requirement that the uncertainty in the final value 
of the energy of motion of H + relative to H must be less than D , 
the binding energy of the final reaction product 

One of the conditions for the use of a classical description 
specifies that the wavelength X of any particle involved must be much 
less than the smallest distance involved in the collision which, in this 
case, is the distance of closest approach r^ , or 

X « r (5) 

a 



For elastic scattering in a central force field, an analysis 

such as that given by McDaniel (1964b) shows that r may be derived 

a 

from the scattering parameter b in the following manner. In the 
center of mass system b is given by the expression 



b 



Z 

M* 



2 

z e 
v 2 



r 



2 



( 6 ) 
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where Z and z are the atomic numbers of the colliding particles, 
is the reduced mass of the system, is the relative velocity 

between the particles, and (g) is the scattering angle. 

For the H + - H collision, 

% 

M 7 and =“ tt/2 (7) 

and for the H - C collision. 



M- “ m H and <g)£ =“ tt/2 . (8) 

Making these substitutions in (6) we obtain 

e 2 6 e 2 

b H+-H = ~ and b H-C = ~E (9) 

where E is the kinetic energy of the incident H + in the laboratory 
frame. Because of different velocities, E is also the kinetic energy 
of the particle after the interaction. 



The distance of closest approach r can be found as the largest 

a 

real root of the equation 



V(r a ) 




( 10 ) 



where V(r ) is the Coulomb potential evaluated at r . Substituting 
a a 

(9) into (10) we obtain 



a H+-H 



2,42 e 2 
E 



and 



r 



a 



H-C 



14,5 e 2 
E 



(ID 
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is the limiting dimension to 



It is therefore seen that r 

a H + -H 

be used for the de Broglie wavelength criterion. To satisfy this 
requirement we take 



\ = — - — = 1/10 r 

"“h V 1 a H+-H 



( 12 ) 



whereby we obtain = 1.69 x 10 cm/sec corresponding to an incident 

kinetic energy of 150 ev. As X ~ and r ^ ~ , it is seen that 

v a 2 

v 

r^ will decrease faster than X , and thus this value represents an 
upper limit on the energy. 

The second requirement for the application of classical theory 
applies the restrictions imposed by the Uncertainty Principle. The 
uncertainty in energy AE introduced by specifying the transverse 
position of the particle must be small compared to the binding energy 
of the final product. If the maximum uncertainty in the transverse 
direction is Ax , and in momentum is Ap , then this requirement can 
be written as 



However as 



we obtain 




AE = 



8m R (Ax) 



(13) 



(14) 



(15) 



and since 



Ax < r 

a 



(16) 
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we have 



AE > 




(17) 



In addition AE should be small compared to D , and we can finally 
write this requirement as 



D » 



8m H 




3 H+-H 



(18) 



using r as the smallest distance of closest approach involved. 

a H + -H 



Again using a factor of 10, we require that 



D 



8m H 



10ft 2 

2 



a 



H + -H 



(19) 



Substituting (11) into (19) and solving for the incident 
energy, we find 

2.42 e 2 V8nU> ' 

E ‘ 



H 



kinetic 



( 20 ) 



where D is 2.6 ev for an particle. Putting in the proper values, 
one calculates E = 820 ev as satisfying this condition. As in the 
de Broglie wavelength criterion, this value also places an upper limit 
on the incident H energy. 



There exists, however, one further peculiarity which should be 
mentioned. Bohm (1951) shows that for a strictly Coulomb potential, 
the classical (Rutherford) cross section is the same as the exact 
Quantum Mechanical cross section even in the region where classical 
approximations do not hold. Therefore although the de Broglie wave- 
length criterion would limit the experiment to be performed for an 
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incident energy below 150 ev, this restriction can be ignored providing 
the scattering events are coulombic, and thus we assume 800 ev as an 
upper limit. 

For a lower limit the restriction is placed that the incident 
particle energy should be greater than any of the binding energies of 
the systems involved. Thus a lower energy limit of about 25 ev should 
constitute a suitable value. 

Therefore, between the limits of 25 and 800 ev there exists a 
region wherein the method of the impulse approximation using classical 
theory is valid. 

(b) Wire Orbit Method 

This method (1960, 1956a) is based on the equivalence 
relating the path of a flexible wire carrying an electric current in a 
magnetic field, and the trajectory of a charged particle passing through 
the same field. 

Consider a small length ds of flexible wire carrying a current 
and under a tension in a magnetic field as shown in Fig. 2. Using 
Gaussian units, we let 

Ti = unit vector normal to current carrying wire 

■jj = curvature of wire (directed away from the center) 

t = unit vector tangent to wire and in direction of current 

i & = current in abamperes 

T = tension of wire in dynes 

m “g = weight per unit length of wire in dynes/cm 

w 

B = magnetic flux density vector in gauss 
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i ds (t x B) 




Fig, 2 



forces on a Small Segment of a Current-Carrying Wire 
in a Magnetic Field 
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It can be seen that ~t , n , and T are coplanar, and at equilibrium 
the sum of the magnetic force, the gravitational force, and the tension 
must equal zero, giving 



i^ds (t x B) - 2T ~n - m^gds = 0 



Letting ds = pda , we have 



i ds (t x *B) - T “ n - m gds = 0 
a p w to 



( 21 ) 

( 22 ) 

(23) 



from 'which the curvature of the wire is obtained in the following form 



n 

P 



(txt) 




(24) 



A similar expression may be obtained for the path of a charged 
particle in the same field. We let 

*” * curvature of particle trajectory 

— 

t = unit vector tangent to trajectory and in direction of velocity 

0 

— = charge to momentum ratio of particle in gaussian units 

The differential equation for the path of the particle is described by 



di = 

dt 



e^ 

c 



(v x B) . 



Now 





ds 



vdt 



therefore 





ev 

c 



(txt) 



(25) 



(26) 



(27) 
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but 



dt _ n 
ds p 

and the curvature of the path is 



( 28 ) 



n 

P 



77 (t x B) . 



(29) 



If the weight of the wire can be neglected , then the curvature 
equations (24) and (29) are identical when 



T 

i 




(30) 



a 

and the wire and the particle path will coincide for equal boundary 
conditions (end points). 



It is convenient to express equation (30) in experimentally 
useable units. We take the charged particle as the scattered , 

e as a unit positive charge, and convert the momentum p to its 
energy equivalent by the expression 



E = 



2 

£_ 

2m 



(31) 



Because of different velocities, the energy of the incoming proton is 
equal to the scattered particle. Therefore, using the appropriate 

conversion units, we have for the expression relating the wire and 
particle paths 

^ . 48 a 

1 



1. A rough calculation made on an average trajectory shows that the ratio 
of the magnetic to the gravitational force is of the order of about 100. 
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where E is the energy of incoming protons in units of electron volts, 
M is the tension of the wire in grams, and i the wire current in 
amperes . 
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3. Magnetic Field Measurements 
(a) Experimental Equipment 

1. Magnet Coil 

The magnet coil (1956c) consists of a set of five concentric 
brass spools wound with number 14 copper wire* For adequate cooling 
1/4 inch o.d. copper tubing is placed at the outer layer of each spool. 
The individual cooling coils are connected to a common manifold and 
then to a cold water supply. Fig. 3 shows the electrical circuit of the 
five layers arranged to distribute the resistances as evenly as possible. 
A rough estimate of the maximum safe coil current indicated a value in * 
the neighborhood of 15-16 amperes. 

2. Power Supply 

In order to fulfill the specifications of an output current 
of 0-15 amperes and a regulation of about + O.l7o to include both line 
fluctuation and voltage ripple, a Sorensen Nobatron model DCR 150-15 
power supply was acquired and found to satisfactory. 

3. Hall Axial Probe 

As the axial component of the magnetic field predominates in 
the area in which we are interested, it was sufficient only to measure 
this component. The principal problem faced in the measurement of the 
field was that of accurately positioning the measuring equipment within 
the long aluminum tube concentric with the axis of the magnetic coil. 
Conventional gaussmeters or nuclear resonance probes could not be 
inserted into the tube at the desired positions. To overcome the above 
problem, a Halltron HP-315 Axial Field Hall probe^ was used. The probe 

1. Ohio Semitronics, Inc., 1205 Chesapeake Avenue, Columbus, Ohio. 
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AA/V 
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18.36 Q 

AA/V 

Coi 1 2 



18.04 0 

AAAr 

Coil 4 



16.80 Q 

AA/V 

Coil 1 



20.86 0 

AW 

Coil 3 



-* T = 10.52 Q 



Fig. 3 

Circuitry of Magnetic Coil 

(Coils 1 through 5 are labeled from inner to outer layer respectively. 
Resistances are measured in ohms.) 
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itself was placed in a plastic capsule with a long 4-wire cable attached 
to the input and output terminals of the probe. Attached to the 
forward tip of the probe was the positioning device consisting of a 
wooden rule accurate to within 1 mm. 

To provide accurate positioning, a 3/4 inch aluminum pipe was 
fitted through two templates as shown in Fig. 4. Cross hairs were 
scribed onto the templates to allow for azimuthal alignment with a 
carpenter's level. Seven 3/4 inch holes were cut, 1/2 inch apart on 
an axis radially from center, to provide for the positioning of seven 
radial measurements for each azimuthal setting of the templates. 

The encapsulated probe was calibrated in a variable magnetic 
field against a known nuclear magnetic resonance (NMR) probe accurate 
at 1000 gauss to +1 gauss. Magnetic field values were calculated for 
proton resonance in water from the expression: 

TT/t . „ \ v (resonant frequency in mc/sec) 

H(kilogauss) = - 2 - 5 y 7 M ^ (32) 

The circuitry consisted basically of steady 100 milliampere current 
fed into the probe with probe output voltages being read across a four 
ohm resistor with a potentiometer readable to +.05 millivolts. Nine 
calibration points were fed into a "least square" curve fitting program 
on the CDC 1604 computer, and a first degree polynomial was found to 
give the best fit. 

(b) Measurements and Results 

The probe system was assembled and mounted on the aluminum 
tube. The initial run consisted of measurements taken 4 cm apart along 
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7 holes 13/16" dia 1/2 
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sltioning of Hall Probt 



the axis of the magnet at current settings of 2, 4, 6, 8, 10, 12, 
and 14 amperes respectively. The complete tabulation of these values 
is to be found in Appendix 1, Table I. A plot of the field strength, 
taken for the above current settings and at fixed locations (z = 0, 14, 42) 
along the magnet axis, is shown in Fig. 5. It can be seen that the field 
strength is linear over the complete range of magnet current settings. 
Plotting the field strength along the axis of the magnet for four 
different current settings (Fig. 6) indicates the symmetry of the field 
along the axis and shows its largest values located at or near the 
z = 0 plane (the center of the magnet). 

In order to determine the axial component of the magnetic field at 
points off the axis, another series of measurements were made in the 
vertical planes defined by z = 0 , z = +14 , z = +42 cm at 1/2 inch 
radial intervals, and at every 45 degrees from 0 to 360. The results 
of these measurements are tabulated in Appendix 1, Table II, and 
plotted in Figs. 1-1 through 1-5. It can be seen that for a given 
radial dimension, the field is azimuthally relatively symmetric. At 
the center of the magnet its strength increases radially approximately 
10 percent, at 14 cm from the center it decreases radially some 7-8 
percent, and at the ends it is constant. 

Upon completion of the magnetic field measurements, the aluminum 
tube was removed as it was not needed for the trajectory work. 
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Fig* 5 

Axial Component of Magnetic Field Strength as Function of Magnet Current 

at Vai‘-\.:s •_ i ons along Ma»»- * 



See Fig., 4 for int ov;r. 
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Axial Component or Magnetic field Strength (gauss) 




Fig, 6 

Axial Component of Magnetic Field Strength along Axis 
for Various Magnet Current Settings 



u e See Fig, 4 for interpretation of z 
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4. 



Particle Trajectory Measurements 



(a) Experimental Equipment and Arrangement 

1. Establishing Magnet Axis 

The use of symmetry suggested that the scattered particles 
of the proposed experiment should both leave and refocus on the magnet 
axis. This dictated that all the measuring equipment used in the wire 
orbit method should be positioned on this axis. Four equidistant 
indentations were scratched into the flange of each side of the magnet. 
Thin strings were stretched across opposite indentations establishing 
two known points along the magnet axis, A surveying transit was placed 
in line with these points some 230 cm distant from the magnet. This 
instrument was selected over a string line or other optical devices as 
it is located outside the orbit area and does not block any of the 
measuring equipment. Rechecking the transit before each run insured that 
all measuring equipment was on line. 

2. Current Carrying Wire 

Drawing from the experience gained at the Lawrence Radiation 
Laboratory of the University of California (1956b), an uninsulated wire^ 
consisting of 14 strands of #44 copper wire was used. It weighed approx- 
imately .08 grams/ft which produced some sag in the trajectory, but not 
enough to cause an appreciable error in the measurements. It was found 
that it could safely take any current up to 4.5 amperes, but that at 
higher currents it would become oxidized and less flexible. 

A few strands of this wire also produced an excellent trajectory, 
but because of the low currents and tensions involved, this wire was 
not used, 

1. Purchased from New England Electric Works, Inc., 367 Main Street, 
Lisbon, New Hampshire. 
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3. Angle Measuring Instrument 

To measure the departure angles of the wire, the measuring 
instrument shown in Fig. 7 was constructed. A copper sleeve was made 
to turn on a one inch o.d. plastic pipe. The sleeve protruded over the 
pipe forming an angle with the pipe axis on which angles from 35 to 55 
degrees were scribed in one degree increments in an arrangement similar 
to that of a micrometer. Those readings which fell between the scribed 
angles were found through interpolation. At the end of the plastic pipe 
a 3/32" o.d. teflon ball was placed on axis and designed to rotate in any 
direction. The wire passed through the pipe and through a small hole in 
the ball and was secured to the rear of the ball by a simple knot. A 
galvanometer was hooked up to the sleeve and indicated contact whenever 
the sleeve touched the current carrying wire. The pipe with its attached 
sleeve was centered inside the core of the magnet along the axis as 
shown in the schematic of Fig. 8. 

A run of fifteen measurements taken for the same angle indicated 
a value accurate to about +30 minutes. Due to the curvature of the wire 
leaving the device, the particle departure angle as measured 0^ does 
not exactly reflect the correct particle departure angle 0^ . An 
expression relating these two angles is derived in Appendix 2 in the 
form of the following transcendental equation: 

R f esc 0 = sin ( f R cot 6 sec 0 ) (33) 

c m c 

where R is the radius of the copper conductor sleeve and f is a 
function of the magnetic field strength and the particle momentum. 

A computer program was written for this expression, with the results 
indicating an additive correction of 30 minutes or less. 
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Positioning of Angle and Tension Measuring Instruments 



4. Tension Measuring Instrument 



One of the objectives of this experiment was to find 



d9 

dB 



the rate of change of particle departure angle with respect to 
the magnetic field at constant energy. A pulley was the immediate 
answer as any weight hanging over it would have kept the ratio of 
M/i constant in the expression 



= 48 t 
1 



independent of manipulation of the magnetic field. Much effort was 
expended into devising a suitable pulley, all of which met with failure 
due to the fact that the inherent friction of the pulley could neither 
be eliminated nor reduced. 

As a result of not finding an adequate pulley, the instrument shown 
in Fig. 9 was constructed to measure the tension in the wire. It con- 
sisted basically of a 4 inch brass-bronze wire shaped in the form o^ an 
"L 11 which followed the principle of a leaf spring. One end was held 
rigid, and the other end contained the flexible wire cia&nped to it by 
a small teflon wedge. The bar of the "L" was placed in contact with 
the bottom of a pointer of an ammeter scale cut down for this purpose. 

Any deflection of the wire spring caused a corresponding deflection on 
the scale pointer. 

This instrument was calibrated by turning it on its side and 
hanging known weights from it. A further check of the calibration was 
made by placing the instrument in its normal upright position and using 
a pulley with the same weights. The calibration curve is shown in Fig. 10. 
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Calibration Plot o£ Tension Measuring Instrument 
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It must be realized that this method of measuring the tension had 



two distinct disadvantages over that of the pulley. First and most 
important of all, there was no way of keeping the tension, and thus the 

energy, constant as one changed the magnetic field. This meant that 

‘ d0 



dB 



much data would have to be taken and the partial derivative of 
could even then only be found indirectly. In addition, as the wire 
spring was displaced, the location where the flexible wire focussed 
back on the magnet axis also moved. This was not too serious as the 
angle of acceptance for the detector to be eventually placed at this 
position is more than large enough to cover this variance. 



5, Orbit Measuring Instrument 

To measure a point on the orbit of the wire we were interested 
in two dimensions; its axial distance z from where the wire left the 
angle measuring instrument, and its radial distance r from the axis, 
as shown in Fig. 8. A simple but extremely functional device (Fig, 11) 
was built to do this. A metric tape was placed parallel to but below 
the magnet axis and referenced to the center of the magnet. A slide with 
an attached pointer was located on a long aluminum angle bar. Supported 
above the slide was an arm, containing the radial measuring device on 
one end, exactly positioned on the magnet axis. A screw on the radial 
measuring device controlled a brass ruler calibrated in mm which could 
be extended until it just touched the flexible wire. The parts of this 
device as well as all other measuring instruments were constructed of 
non-magnetic materials. 
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6. Positioning of Measuring Equipment 

It was experimentally determined that the best wire trajectory 



occurred if the angle measuring instrument was placed on the magnet axis 
with the wire departing near the center of the magnet where the field 
has its greatest strength. In this position the wire could depart at 
any angle between 35 and 55 degrees. The tension measuring device was 
positioned approximately 40 cm down axis from this point as shown in 
Fig. 8. A closer positioning would cause the wire to touch the sides 
of the magnet for large departure angles, and for a more distant position, 
the wire would have too much sag. Fig. 12 shows the wire trajectory and 
all measuring equipment in place. 

(b) Results and Discussion 

1. Existence of a Unique Trajectory 

Consider a charged particle leaving a point in a uniform 
magnetic field with its initial velocity vector at an angle 9^ with 
respect to B . This particle will describe a helical trajectory and 
return to the original field line a distance L from its starting point. 
It is shown in Appendix 3 that the following relationship holds. 



LB = KE^ cos 9 



(34) 



c 



where L 



axial distance travelled in one revolution 



B 



magnetic field strength 



K 



constant 



E 



particle energy 



9 



departure angle 



c 
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Equation (34) shows that for given values of L , E , B and 9^ there 
exists a unique trajectory, since both the force on the particle and the 
initial conditions are completely fixed by (34)* 

Since it is not possible to write down the corresponding expression 
for the case of non-uniform fields, a qualitative experimental test was 
carried out as follows: 

Keeping the field current and the length along the axis between 
measuring instruments constant, but varying the actual length of wire, 
it was noticed that both the equivalent particle energy and the 
departure angle changed. This experiment was performed several times, 
and with each trial a different trajectory was produced. It seems then 
that there is indeed a unique trajectory when the parameters 9^ , B , 

L , and E are specified. 

2. Particle Trajectory Measurements 

It was originally hoped that in the collision experiment one 
could analyze the angular distribution of the scattered particles 
leaving at a given energy by using baffles and by adjusting the magnetic 
field to examine a narrow angular spread at a particular angle. To 
determine the feasibility of this idea, the wire orbit method was utilized. 

The first objective was to examine the actual particle trajectories 
for different values of departure angle, particle energy, and magnet 
current. A complete tabulation of this data is to be found in Appendix 4. 
In order to have some feeling for the general shape of these trajectories, 
two plots have been made. Fig. 13 shows two trajectories for approx- 
imately equal particle energies but at different angles, while Fig. 14 
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Fig 13 

Plot of Particle Trajectories for Approximately Equal Particle Energies 
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Fig. 14 

Plot of Particle Trajectories for Approximately Equal Departure Angles 



shows three trajectories for approximately equal angles but at different 
particle energies. Although the actual trajectory follows the path of 
a spiral, only the radial distance of the trajectory from the magnet 
axis has been plotted. Note that each trajectory starts off linearly, 
reaches a maximum curvature between 9 and 14 cm from its initial 
departure point, and then returns linearly to the axis. 

In the second objective, we were interested in determining the 

rate of change of particle departure angle with respect to the magnet 

current while holding the particle energy constant. This data is 

tabulated in Table I and plotted in Fig. 15. Because this data had to 

be obtained indirectly, it is somewhat less reliable than the directly 

measured quantities. However, in the 200-500 ev range the value of 
f \ 

— seems quite consistent and has a value of about 3 /amp. 
i B/E 
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Table I 



Calculation of Rate of Change of Particle Departure Angle 
with Respect to Magnet Current at Constant Energy 



Corrected angle a 
of departure 
(degrees) 


Particle 

Energy 

(ev) 


Magnet 

Current 

(amp) (d< 


fd 0 | 

KJe 

sgrees/amp) 


52.74 


108 


7 


10.09 


42.65 


108 


8 




48.12 


133 


8 


9.37 


38.75 


137 


9 




45.44 


193 


9 


3.59 


41.85 


193 


10 




45.80 


338 


11 


2.68 


43.12 


332 


12 




45.95 


440 


12 


2.85 


43.10 


445 


13 




47.09 


630 


13 


0.76 


46.52 


640 


13.75 





a. All angles are accurate to + 0.50 degrees. 
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Rate of Change of Particle Departure Angle with 
Respect to Magnet Current at Constant Energy 



5. 



Cone lusion 



Through the wire orbit method it has been shown that by using the 
existing magnet, it is indeed possible to obtain trajectories in 

the required energy range of 25-800 ev having departure angles between 
35 and 55 degrees. For the above energy and angular requirements, it 
was experimentally found that one set of trajectories exists which 
depart near the center of the magnet and return to the axis approximately 
40 cm away. Although valid trajectories may exist for other end point 
positions, time did not permit these to be investigated. 

From Figs. 13 and 14, one concludes that for fixed initial and 
end points, the trajectories are very similar to each other. This 
would seem to make the use of baffles as a particle trajectory 
discriminating mechanism very doubtful. 



Typical values of the rate of change of departure angle with respect 

'd0 
dI T 



B 



were found to be of the 



to magnet current at constant energy 

order of 3 deg/amp in the energy range of 200-500 ev. 

It is expected then that the technique to be utilized for 
differentiating particles emerging with different departure angles at 
a given energy will involve both moving the detector along the axis 
and varying the magnet current. This paper has shown that both these 
methods can be used and the numerical analysis presently being conducted 
should produce the necessary quantitative information needed. 

Finally, it is shown that, if the particle energy, departure angle, 
magnetic field, and the distance between the initial and end points of 
the trajectory are fixed, there exists a unique trajectory. 
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Appendix I 
Table I 



Axial Magnetic Field Strength along Axis of Magnet for Various Currents 



Distance along Magnetic Field Strength (gauss) 

Axis from Center 



of Magnet (cm) 


1=2 


4 


6 


8 


10 


12 


14 amp 


z a = 66 


7 


12 


16 


18 


23 


28 


30 


62 


9 


14 


18 


21 


25 


30 


37 


58 


10 


16 


20 


25 


31 


38 


42 


54 


10 


18 


23 


31 


36 


43 


47 


50 


12 


20 


29 


34 


43 


49 


57 


46 


14 


23 


32 


43 


53 


59 


69 


42 


16 


29 


40 


53 


64 


72 


85 


38 


20 


34 


49 


65 


80 


92 


105 


34 


25 


43 


62 


84 


100 


119 


137 


30 


31 


56 


82 


108 


131 


152 


176 


26 


40 


73 


106 


141 


172 


203 


234 


22 


51 


95 


142 


186 


232 


269 


313 


18 


67 


128 


190 


254 


315 


373 


430 


14 


89 


172 


258 


342 


429 


502 


590 


10 


117 


229 


342 


454 


568 


675 


795 


6 


144 


284 


425 


564 


704 


843 


985 


2 


163 


319 


478 


634 


793 


945 


1110 


0 


164 


322 


484 


641 


805 


958 


1125 


- 2 


163 


319 


478 


632 


793 


945 


1110 


- 6 


144 


284 


421 


562 


700 


843 


985 


-10 


117 


227 


341 


454 


564 


675 


795 


-14 


89 


172 


256 


341 


427 


502 


590 


-18 


67 


128 


190 


252 


311 


373 


430 


-22 


51 


93 


141 


186 


230 


269 


313 


-26 


38 


71 


104 


139 


172 


203 


234 


-30 


29 


54 


82 


108 


131 


152 


176 


-34 


23 


43 


62 


84 


100 


119 


137 


-38 


18 


34 


49 


64 


80 


92 


105 


-42 


16 


27 


40 


53 


64 


72 


84 


-46 


14 


23 


32 


43 


53 


59 


69 


-50 


12 


20 


27 


36 


43 


49 


57 


-54 


10 


18 


25 


31 


36 


43 


47 


-58 


10 


16 


20 


27 


31 


38 


42 


-62 


9 


14 


18 


21 


27 


30 


37 


-66 (Detector End) 


7 


12 


16 


18 


25 


28 


30 



a. See Fig. 4 for orientation of coordinate system. 
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Appendix I 
Table II C 

Axial Magnetic Field Strength off Axis of Magnet 



Distance Radial 

along Axis Distance Magnetic Field Strength (gauss) 

from Center from Magnet 



of Magnet (cm) 


Axis (cm) 


o 

II 

Q 

CD 


45 


90 


135 


180 


225 


270 


315 


Z b = 42 


r b = 0.00 


62 


64 


64 


62 


62 


64 


64 


62 




1.27 


62 


64 


64 


62 


62 


62 


64 


62 




2.54 


62 


64 


62 


62 


62 


64 


62 


62 




3.81 


62 


62 


62 


62 


62 


64 


62 


62 




5.08 


62 


62 


62 


62 


62 


62 


62 


62 




6.35 


62 


62 


62 


62 


60 


62 


62 


62 




7.62 


62 


62 


62 


62 


60 


62 


62 


60 


z = 14 


r = 0.00 


427 


429 


429 


429 


429 


429 


429 


429 




1.27 


425 


427 


427 


427 


429 


429 


429 


429 




2.54 


421 


421 


425 


425 


427 


429 


427 


427 




3.81 


418 


418 


419 


419 


427 


423 


423 


419 




5.08 


412 


412 


412 


414 


421 


421 


419 


410 




6.35 


401 


399 


401 


407 


410 


421 


414 


407 




7.62 


396 


385 


388 


394 


399 


405 


403 


399 


N 

II 

O 


r =0.00 


801 


801 


805 


803 


803 


805 


803 


803 




1.27 


801 


806 


806 


806 


805 


805 


805 


805 




2.54 


806 


808 


806 


810 


810 


814 


806 


814 




3.81 


823 


823 


823 


823 


819 


823 


823 


823 




5.08 


843 


843 


843 


843 


841 


843 


838 


843 




6.35 


861 


863 


860 


867 


860 


861 


861 


865 




7.62 


889 


891 


887 


893 


887 


894 


887 


896 


z = -14 


r =0.00 


427 


425 


427 


427 


427 


427 


425 


427 




1.27 


427 


427 


425 


425 


421 


427 


427 


427 




2.54 


425 


427 


421 


421 


418 


423 


419 


423 




3.81 


419 


423 


419 


416 


414 


419 


418 


419 




5.08 


414 


418 


414 


410 


405 


410 


408 


410 




6.35 


405 


412 


405 


403 


396 


401 


401 


405 




7 .62 


394 


403 


397 


394 


385 


388 


388 


390 


z = -42 


r =0.00 


64 


64 


64 


62 


62 


64 


64 


62 




1.27 


62 


64 


64 


62 


62 


62 


64 


62 




2.54 


64 


64 


64 


62 


62 


64 


64 


62 




3.81 


62 


62 


62 


62 


62 


64 


62 


62 




5.08 


62 


62 


62 


62 


62 


62 


62 


62 




6.35 


62 


62 


62 


62 


60 


62 


62 


62 




7.62 


62 


62 


62 


62 


60 


62 


62 


60 



a. All values of magnetic field strength were measured at a constant 
magnet current of 10 amperes. 

b. See Fig. 4 for orientation of coordinate system. 

c. Figs. 1-1 through 1-5 show this data plotted. 
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Fig o 1-jl 
( z = +42 cm ) 
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Fig, 1-2 
(z - +14 cm) 
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90 ° 




Fig. 1-3 
(z = 0) 
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90 ° 




F-g 

(z = - 14 cm) 
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Fig, 1-5 
(z = 
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Appendix 2 



Derivation of Expression Relating Corrected and Measured 
Particle Departure Angles 



The charged particle leaves the magnet axis as shown in Fig. 2-1 
with an initial angle 9^ and velocity v . Assuming that near the 
center of the magnet where the particle leaves, the magnetic field 
strength is uniform, then the particle assumes a radius of curvature 
p as given by 



P 



m v c sin 6 

c 



e B 



(2-1) 



where all units are gaussian, and m is the particle mass, v is 
the initial particle velocity, 9^ is the actual particle departure 
angle, e is the charge of particle, B is the magnetic field strength 
in gauss, and c is a dimensionless factor 3 x 10^. Because its 
velocity is skew to the magnetic field, the particle traverses a spiral 
whose radial coordinate r in plane polar coordinates is 



r = 2 p sin cp 



( 2 - 2 ) 



and whose axial component is given by 



z = k cp 

where k is the pitch of the spiral. 

To find k we know that in one period T where 

T _ 2 TT m c 

e B 



(2-3) 



(2-4) 
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Schematic Showing Wire Departure from Angle Measuring Device 



the particle passes through cp = tt radians and travels an axial 
distance equal to 



z = v cos 9 T 
c 



(2-5) 



Putting (2-4) into (2-5) and equating this to (2-3) with cp = TT , 
we find that the pitch k is expressed as 



k = 



2 m v c 
e B 



cos 



( 2 - 6 ) 



Using (2-6) we solve (2-3) for cp in terms of z with the 
result that 



cp = 



z e B 

2 m v c cos 9 



(2-7) 



Substituting (2-1) and (2-7) into (2-2) we have 



r = 



2 m v c 
e B 



sin 9 sin 



e B z 

2 m v c cos 9 



( 2 - 8 ) 



Applying the boundary condition that when the wire touches the 
conductor of the angle measuring device 



r = R and z = R cot 0 



m 



(2-9) 



where 9^ is the measured particle departure angle, 



Letting 



f = 



e B 



2 m v c 



(2-10) 



we rearrange terms and derive' the following transcendental equation 
relating the corrected and measured particle departure angles as 



R f esc 0 = sin (f R cot 0 sec 0 ) . (2-11) 

c me 
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Appendix 3 

Conditions for Unique Trajectory in Uniform Field 



From the differential equation describing the motion of a charged 
particle in a uniform magnetic field it is easy to show that 

= e vj. B (3-1) 

where m , e , P , and vj_ is the mass, charge, radius of curvature, 
perpendicular component of the velocity of the charged particle 
respectively, and B is the magnetic field strength. Eliminating vx 
from the right side of (3-1) and writing 

m vx = m v sin 9^ (3-2) 

where v is the magnitude of the velocity of the charged particle 
and 9^ the initial angle that it makes with the field, we have 

v sin 9 = e - - - P - . (3-3) 

cm 



The period T of the charged particle is given by the expression 



T = 2 17 P 

v sin 9 

c 

from which, upon substituting (3-3) into (3-4), we get 

_ 2 tt m 

1 e B 



(3-4) 



(3-5) 
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For a given length L along the axis of the field, the charged particle 
will refocus back on the axis in one period providing 



Substituting 



L = v cos 9 T 
c 



}-5) into (3-6), we have 



2 tt m v 



cos 9 



(3-6) 



(3-7) 



Writing this expression in terms of the particle energy, we have 



L B 



9 3/2 tt * 
2 TT m 



cos 



(3-8) 



As the right hand side is a constant, we can write (3-8) as 

L B = K cos 6 (3-9) 

c 

which defines any one of the four parameters L , B , E , and 0^ 
in terms of the remaining three. This then is the condition that 
determines a unique trajectory over one period in a uniform magnetic 
field. 
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Appendix 4 

Particle Trajectory Measurements from Wire Orbit Method 



Here 6 is the measured departure angle, 0 is the corrected departure 
m c 

angle, m is the tension, i is the wire current, E is the energy 
of the particle, I is the magnet current, z is the distance along 

the magnet axis, and r is the radial distance from the magnet axis. 



0 a 


9 J 


m 


i 


E 


I 


z a 


r a 


m 


c 














(degrees) 


(degrees ) 


(gms) 


(amp) 


(ev) 


(amp) 


(cm) 


(mm) 


45.48 


45.78 


0.86 


2.0 


75 


6 


2.5 


26.0 














4.5 


51.0 














6.5 


69.0 














8.5 


77.0 














10.0 


79.0 














12.5 


79.0 














14.5 


78.0 














17.5 


74.0 














20.5 


66.0 














24.5 


55.5 














27 .5 


45.5 














30.5 


34.5 














33.5 


24.5 














36.5 


12.5 














40.4 


0.0 


47.00 


47.25 


0.48 


2.5 


84 


6 


2.5 


27.5 














4.5 


52.5 














6.5 


70.5 














8.5 


77.5 














10.5 


80.0 














12.5 


79.5 














14.5 


77.5 














16.5 


74.5 














19.5 


67.0 














23.5 


56.5 














26.5 


47.0 














29.5 


36.5 














32.5 


26.0 














35.5 


14.5 














40.4 


0.0 



a. See Fig. 8 for orientation of coordinate system. 
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0 

m 


e 

c 


m 


i 


E 


I 


z 
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47.85 


48.10 


0.57 


3.0 
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2.5 
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